Take-off of small Leidenfrost droplets 
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We put in evidence the unexpected behaviour of Leidenfrost droplets at the later stage of their 
evaporation. We predict and observe that, below a critical size Ri, the droplets spontaneously take- 
off due to the breakdown of the lubrication regime. We establish the theoretical relation between 
the droplet radius and its elevation. We predict that the vapour layer thickness increases when the 
droplets become smaller. A satisfactory agreement is found between the model and the experimental 
results performed on droplets of water and of ethanol. 



PACS numbers: 47.55D-,68.03.-g 

The name of J.G. Leidenfrost (1715-1794) is still re- 
membered because he was the first to publish the ob- 
servation that a puddle of water dropped on a very hot 
surface divides into droplets, which stand each above the 
surface and slowly evaporate. At the end of the pro- 
cess, he noticed [1] : "... before the whole drop dis- 
appears. Which at last exceedingly diminished so that 
it can hardly any more be seen, with an audible crack, 
which with the ears one easily hears, it finishes its exis- 
tence, and in the spoon [making the hot surface] it leaves 
a small particle of earth ..." 

Over the years, the Leidenfrost phenomenon has at- 
tracted the attention of many investigators, see for in- 
stance 0-05 seemingly little interest has been paid 
to the observation by Leidenfrost of the final stage of the 
droplet evaporation, in particular what happens when 
the droplets become very small before finally disappear- 
ing. We investigate in this letter theoretically and ex- 
perimentally the final regime of the Leidenfrost droplets. 
We show that when the droplets become smaller than 
a well defined radius, they suddenly take-off due to the 
breakdown of the classical Leidenfrost lubrication regime. 
They thus reach an elevation which is much higher than 
their radius. 

Let us consider first an almost spherical drop of volatile 
liquid with radius R and mass M standing slightly above 
a hot flat surface with a vapour film of thickness h in 
between. Let I be the horizontal width of this film, this 
length will be related to the other physical parameters of 
the problem. 

Our estimate will rely on the fact that, in a range of 
parameters, the horizontal length scale I is much larger 
than h. Both I and h are less than R, the radius of 
the droplet. Furthermore, we shall consider a range of 
parameters where the vapour pressure in the film is neg- 
ligible compared to Laplace's capillary pressure. There- 
fore the droplet remains almost exactly spherical in this 
regime. This will yield the condition (j3]) derived below, 
valid under the constraint h <C R. This condition hap- 
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FIG. 1: Experimental set-up. A cloud of sub-millimetric 
droplets is rained toward a silicon wafer inserted in a cop- 
per bloc kept at a controlled temperature Ti. The droplet 
of radius R is at temperature To and at a distance h of the 
substrate. 



pens to be impossible to satisfy if the droplet becomes 
too small, that is if its weight is too small to balance the 
upward pressure in the gap between the hot plate and 
the droplet. Therefore, it is natural to guess that, when 
this happens, the droplet lifts off the plate. This defines 
a second regime where h ^> R. The scaling laws for this 
second regime are given in this letter and we show that 
the droplet reaches an higher elevation as it gets smaller 
and smaller by evaporating. 

Energy conservation during the evaporation process 
of the droplet (namely Stefan's boundary condition on 
the liquid/vapour interface), which we assume to be just 
below the boiling temperature, and Fourier law for the 
heat transfer in the film yield the order of magnitude for 
the vertical velocity of vapour w near the surface of the 
droplet: 



w = 



X5T 
hLp v 



(1) 



Here 5T = T\ — To is the temperature difference between 
the droplet and the hot plate, p v is the mass density of 



the vapour, A the heat conductivity in the vapour and 
L the latent heat of evaporation per unit mass. Lubri- 
cation theory, valid if I ^> ft, yields the magnitude of 
the horizontal velocity u in term of the vertical velocity 
: u = where I is the horizontal extent of the vapour 
film. From Poiseuille relation such a flow is driven by 
a horizontal pressure gradient scaling as SP/l ~ r]u/h 2 . 
Therefore the typical gradient is of order SP/l ~ wlr]/h 3 . 
Using now equation (pQ) we find (replacing sign ~ by = 
for readability) : 



SP = r]l 2 



XST 
h 4 Lp v ' 



(2) 



The liquid drop is at mechanical equilibrium when 
Mg = 5 PI 2 namely when its weight is equal to the ver- 
tical component of the force due to viscous pressure in 
the gap between the droplet and the hot plate. Note that 
the contribution of the viscous stress to this vertical force 
is of the same order of magnitude as the one of hydro- 
static pressure. Therefore we find using the mechanical 
equilibrium relation given above that : 



V_ = Lp v Mg 
ft 4 rjSTX ' 



(3) 



As stressed above the lubrication approximation be- 
comes invalid when ft ~ I. From the M dependence of 
the right-hand side of equation (j3j) we find that the lubri- 
cation approximation becomes invalid if the radius of the 
droplet is less than a critical drop radius Ri defined by 
the condition I = ft. The breakdown of the lubrication 
approximation will then occur for : 



Ri 



( rjSTX 

\gLp v pi 
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(4) 



Therefore one expects that, as R gets smaller than R\, 
the droplet takes off from the hot plate. The range of 
values of ft ~ R ~ Ri is beyond reach of simple order of 
magnitude estimates and its analysis requires a solution 
of the full set of equations for the multi-phase Stokes flow 
and the heat transfer model. Let us consider the asymp- 
totic regime where R is much less than R\ , where one ex- 
pects the droplet to lift well above the hot plate (ft ^> R). 
In this regime, the equilibrium elevation of the droplet is 
reached when the gravity force is balanced by the dipolar 
component of the velocity field generated by the evapo- 
ration from the droplet surface and by its image with 
respect to the horizontal plane. This dipolar approxi- 
mation is correct for ft ^> R. Let T\ the temperature 
at infinity and To the temperature of the droplet. The 
temperature field satisfying Laplace's equation around 
an isolated droplet reads : T(r) = T x + (T - T x )R/r. 
The origin of the spherical coordinates (r = 0) is taken 
at the center of the droplet. Let us consider now the 
temperature due to the image with the temperature field 
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FIG. 2: Top : Pictures of a Leidenfrost droplet and its re- 
flected image on the silicon wafer. Bottom : Height of the 
droplet as a function of time. After being dropped, it first 
bounces and lands on the substrate. The droplet then stands 
on the substrate in the lubrication regime. It finally takes-off 
from the substrate. 



of order T^ m ~ (T\ — To)R/2h near the physical droplet, 
with a nearly uniform vertical gradient. As justified be- 
low, in the limit of large ft the net vertical heat flux is of 
order of the derivative of T im with respect to ft, namely 
(Ti — To)R/h 2 . Using Stefan's law we find the order of 
magnitude of the net vertical velocity in the vapour near 
the droplet reads : 



A (Ti-To)iJ 



h 2 



(5) 



This relation is similar to equation (pp), but with ST/h 
replaced by (Ti - T )R/h 2 . Multiplying this velocity by 
r]R, as in Stokes drag law, we obtain the vertical force 
on the droplet. Balancing this force with gravity we find 
that, in the limit h ^> R (namely for R <C i?/), the 
dimensionless elevation of the droplet at equilibrium h* = 
h/R should scale as : 



ft* = (Ri/R) 



3/2 



(6) 
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Before going further we have to analyse the constraint 
on the pressure in the gap: it must be negligible com- 
pared to the Laplace's pressure in the droplet to ensure 
that it remains close to its spherical shape everywhere. It 
can be shown [10] that I ~ V Rh is a fair estimate of the 
horizontal extent of the layer for which the viscous pres- 
sure contributes significantly to the vertical force. Using 
eq.(j2j) and this relation between Z, R and ft, the pres- 
sure induced by the Poiseuille flow in the gap is found 

to be of order Sp = gpi^§r- It crosses the order of 
magnitude of Laplace's pressure inside the droplet when 
SP = j^. This happens when the radius becomes of or- 
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(^f^j - This relation is obtained using 

the above geometrical relation for Z and equation (|3j) re- 
formulated as ft// = (R1/R) 3 / 4 . This assumes that the 
/ \ 1/2 

capillary length ( J is much larger than R^ since 

otherwise the cross over would happen in a range of val- 
ues of R smaller than Ri'. this is impossible because Ri 
was derived under the assumption that the lubrication 
approximation holds in the gap between the sphere and 
the hot plate. The condition Ri ^> Ri is fulfilled in the 
experiments described below, but there could be other 
situations at stakes for which this condition may not be 
fulfilled. This gives us also the opportunity to address an- 
other feature of our experiments, namely the difference 
of behaviour between droplets of water and of ethanol. 
Since the quantity Ri / Ri is independent on the radius of 
the droplet, but depends on the properties of the liquid 
and its vapour, the behaviour of Leidenfrost droplets is 
not unique up to rescalings, because it depends on a di- 
mensionless ratio with different value depending on the 
liquid used. Let us also notice that for radii larger than 
Ri, the shape of the droplet becomes rather complex, be- 
cause it depends on the (large) deformation of the droplet 
surface due to the Poiseuille pressure inside the vapour 
gap [ll|, ll2| . This happens in a range of parameters out- 
side of this study. Therefore we focus on the range of 
parameters such that Ri <C R <C Ri. In this range, us- 
ing equation (|3]) and the geometrical relation for Z, the 
dimensionless height satisfies: 



ft* = (Rt/R) 



3/2 



(7) 



For usual liquids the value of Ri is of a few hundred 
of microns. It is therefore worth mentioning that this 
prediction is qualitatively different and even opposed to 
the one commonly used for which the vapour thickness 
is expected to decrease when the radius of the droplet 
decreases. 

The previous relations rely on the assumption that the 
the temperature field between the drop and the substrate 
(the plane located at z = 0) is a solution of Laplace's 
equation. This approximation is valid in the limit in 
which the convective heat flux is negligible compared to 



the diffusive flux (small Peclet number limit) and it holds 
in our experiments. From the image solution of potential 
problems, the temperature field between the drop and the 
hot plate is the same than the one between two spheres at 
temperatures To and 27\ — To and separated by a distance 
2ft. Laplace's equation for two spheres was solved one 
hundred years ago by Jeffery [13] thanks to bi-spherical 
coordinates. We have numerically recovered the expected 
dependence of the temperature gradient in the limits of 
h/R large and h/R small. As given above by scaling 
arguments, they respectively read VT ~ ATR/h 2 and 
VT-AT/ft. 

Considering a water drop on a hot plate at 400°C the 
numerical estimate of the critical radius is R\ ~ 19/im. 
To our knowledge no quantitative experiment has been 
done at such small sizes. We therefore decided to perform 
experiments on ultra-distilled water and ethanol drops 
to verify the scaling laws derived just above. The ex- 
perimental set-up is depicted in Fig. [TJ A cloud of sub- 
millimetric droplets is rained on a silicon substrate kept 
at a controlled temperature T\. The two upper blocks 
are necessary to reduce the number of drops falling in the 
field of the camera. They are also useful for the thermal 
stability of the system. Experiments are performed un- 
der a laminar air flow protection in order to prevent dust 
contamination. A high-speed camera, at a frame rate of 
2000 frame per second, was used to record several evapo- 
ration take-off processes. An image analysis allowed us to 
measure the radius and elevation of droplets as a function 
of time. The typical error made on these measurements 
is 2/xra. On the top of Fig. [2j we show six pictures taken 
during the take-off of a water droplet. On the bottom 
of Fig. [2] we plot the height of the water droplet de- 
posited on a silicon substrate at T = 375°C and room 
pressure. It first bounces and lands on the silicon sub- 
strate. It then enters in the lubrication regimes in which 
both radius and vapour thickness decrease with time. At 
certain time the droplet takes-off from the substrate (see 
the movie Ml in the supplementary materials). 

We recorded several evaporation processes and we rep- 
resent in Fig. [3]the dimensionless elevation of the droplet 
ft* as a function of its radius R. The measurements are in 
qualitative agreement with the pictures of Fig. 2 (top) 
: almost all droplets take-off from the substrate when 
their radii becomes of order Ri. The theoretical predic- 
tion given above is tested against experiments. The full 
line is a best fit with a single parameter a multiplying the 
predicted critical radius Ri . We can see that the theoret- 
ical prediction fits rather well the experimental data with 
a value a = 2.2 close to unity. One can see that there is an 
important dispersion in sizes at which the droplets take- 
off. It seems that, for the system of water droplets under 
study, the lift-off is analogous to a sub-critical transi- 
tion (in other terms for the same radius and temperature 
difference, there are two equilibrium elevation in a cer- 
tain range of values of R/Ri). The same behaviour has 
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been observed for water droplets on substrates at differ- 
ent temperatures. Decreasing the temperature tends to 
decrease the dispersion of radius's at take-off. Neverthe- 
less the discontinuous character of the escape from the 
lubrication regime remains. 

We investigated also the behaviour of ethanol Leiden- 
frost droplet. The same experimental procedure is used 
and the silicon substrate is kept at similar temperatures 
than for the study of water Leidenfrost droplets. Con- 
trarily to water droplets almost all ethanol droplets do 
not rebound and land on the substrate (see the movie 
M2 in the supplementary materials). This is due to the 
lower latent heat of ethanol and therefore to its larger 
evaporation rate. As a consequence the ethanol droplets 
enter directly in the regime described above for which 
h » I and we do not observe a transition from the lu- 
brication regime. The experimental results are presented 
in Fig. [4] in which we plot, as for water droplets, the di- 
mensionless elevation of the droplet h* as a function of its 
radius R. The data are presented for 10 different ethanol 
droplets on a substrate kept at T = 400° C. As stressed 
just above the escape from the lubrication regime is no 
longer present and the height h of the droplets increases 
in a continuous manner. In the insert of Fig. H]we repre- 
sent the same experimental data but plotted in a log-log 
representation. The black line corresponds to the scaling 
law predicted in equation (j6]). A rather good agreement 
is found between the theoretical prediction and the ex- 
periments. 

In this letter we have put in evidence the unexpected 
behaviour of a Leidenfrost droplet in the latter stage of 
its evaporation. Below a critical size Ri the droplets es- 
cape from the lubrication regime and take-off from the 
substrate. In an intermediary regime (Ri <C R <C Ri) 
the thickness of the vapour film increases as the radius 
of the droplet decreases. This latter prediction, based 
upon a detailed analysis of the various physical phenom- 
ena involved, is qualitatively different, and even opposed 
to the one commonly used: as a result of the standard ap- 
proach to the Leidenfrost phenomenon the vapour thick- 
ness is predicted to decrease as the radius of the droplet 
decreases. The different scaling laws presented in this 
paper are in semi-quantitative agreement with the exper- 
imental data obtained on water and ethanol Leidenfrost 
droplets. In particular, as observed and predicted, the 
thickness of the vapour film increases as the droplet be- 
come smaller. Besides its fundamental interest this study 
should find implications in many domains such as for ex- 
ample in diesel combustion engines or in heat transfer us- 



ing cooling spray [14]. We hope that this study will moti- 
vate numerical approaches complementary to the present 
study. 
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FIG. 3: Dimensionless height of the droplets as a function 
of their radii. The red line is a best fit to the equation ©. 
Note that a minority of droplets fall back on the substrate. 
This is certainly due to dust contamination. The same data 
are represented in the insert but in log-log. The dotted line 
indicates the value of Ri for water Leidenfrost droplets on a 
substrate kept at T = 400° C. 
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FIG. 4: Dimensionless height of ethanol droplets as a function 
of R. The substrate is kept at T = 400° C. The same data 
are represented in the insert but in log-log. The black line 
corresponds to the scaling law predicted in equations ©. 



